#ifndef __GAUSS_JACOBI_H__
#define __GAUSS_JACOBI_H__

#include "matrix.h"
#include "vector.h"

namespace t2c2 {
	
class Gauss_Jacobi {
public:
	Gauss_Jacobi(Matrix*, const Vector&);
	~Gauss_Jacobi();

	Vector solve(const Vector&);
	double sum(const Vector&, unsigned int);
	
private: //[M]*x=b
	Matrix* matrix;
	Vector b;
};

}

#endif
